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We present a novel extension of the concept of scars for 
the wave functions of classically chaotic few-body systems 
of identical particles with rotation and permutation sym¬ 
metry. Generically there exist manifolds in classical phase 
space which are invariant under the action of a common sub¬ 
group of these two symmetries. Such manifolds are associated 
with highly symmetric configurations. If sufficiently stable, 
the quantum motion on such manifolds displays a notable 
enhancement of the revival in the autocorrelation function 
which is not directly associated with individual periodic or¬ 
bits. Rather, it indicates some degree of localization around 
an invariant manifold which has collective characteristics that 
should be experimentally observable. 

PACS numbers: 03.65.Sq, 05.45.-|-b 

During the last decade, spectral fluctuation proper¬ 
ties of classically fully chaotic systems (“quantum chaos” 
for short) have been largely understood, and interest has 
shifted to the non-generic properties of wave functions. 
In systems with few degrees of freedom, these functions 
may display scars m, *.e. possess increased intensities 
along short and not very unstable periodic orbits of the 
corresponding classical system. 

In this Letter, we extend the study of scars to systems 
with identical particles. Such systems possess symmetry 
properties which are generically absent in systems with 
two or three degrees of freedom. Particle identity leads 
to permutational symmetry. In addition, a self-bound 
system possesses rotational symmetry. We show that the 
combination of these two symmetries leads to a novel 
mechanism for the formation of scars. 

In a previous paper, two of the authors studied the in¬ 
fluence of both symmetries on the periodic-orbit struc¬ 
ture of chaotic systems [||. In the present Letter, we 
proceed differently. The combination of rotational and 
permutational symmetry allows us to construct invariant 
manifolds in classical phase space. We show that such 
manifolds may lead to scarred wave functions of the cor¬ 
responding quantum system. The scars are not related 
to individual periodic orbits but are similar to the ones 
found by Prosen in billiards 0. We believe that our re¬ 
sults are important for finite many-body systems such as 
atoms, molecules or atomic nuclei where this generalized 
scarring has obvious implications for collective motion. 
Our results derive from the numerical study of a specific 
Hamiltonian. For reasons given below, we believe our 
results to be generic, however. 

The paper is organized as follows. We define a ro¬ 


tationally invariant few-body system of four interacting 
identical particles in two dimensions. We prove the ex¬ 
istence of an invariant manifold and analyze its classical 
stability. Finally, we investigate the corresponding quan¬ 
tum system. This is done by semiclassical propagation 
of wave packets. 

Our Hamiltonian has the form 

where and = {x^,yi) with i = 

are the two-dimensional momentum and position vec¬ 
tors of the particle, respectively. We use units where 
m = a = l,h = 0.01; then coordinates and momenta 
are given in units of and , 

respectively. For the Hamiltonian, this leads to the scal¬ 
ing relation iL(72p,y4r) = 'yH{p,r). This shows that 
the structure of classical phase space is independent of 
energy. Moreover, energy and total angular momentum 
are the only integrals of motion, and the system is non- 
integrable. 

Construction of the invariant manifold. The Hamilto¬ 
nian ( 0 ) is invariant under the symmetric group S'4 and 
the orthogonal group 0(2). Thus, we may apply the 
ideas of ref. Q . By imposing suitable symmetry require¬ 
ments onto the initial conditions, we can restrict classical 
motion to a low-dimensional invariant manifold. Many 
such manifolds exist and it suffices to consider one of 
them. We choose the initial conditions in such a way 
that positions and momenta exhibit the symmetry of a 
rectangle. Such a configuration is shown in Fig. 0 Po¬ 
sitions and momenta of the particles are indicated by 
points and arrows, respectively. Obviously, particle 3 is 
the image of particle 1 under inversion, whereas particle 
2 and particle 4 are the mirror images of particle 1 un¬ 
der a reflection at appropriately chosen axes, respectively. 
The manifold is spanned by the two two-dimensional vec¬ 
tors p = (j>x,Py) and r = {x,y) giving the momentum 
and position of particle f. The associated Hamiltonian 
H{p,r) = ip^ -I- has been studied extensively in 

the literature, both in the classical and the quantum case 
[p[-p^. The classical system is essentially chaotic. Only 
one stable periodic orbit is known which is surrounded by 
a very small island of stability ||^. The periodic orbits 
can be enumerated by means of an (incomplete) symbolic 
code consisting of three letters |]^ . Application of rota¬ 
tions and/or permutations to the configuration shown in 
Fig. 0 yields further equivalent invariant manifolds. 


f 


Stability of motion close to the invariant manifold. For 
particle positions far from the origin, we may expand the 
potential of the Hamiltonian (H) around the “channel” 
configuration defined by the position x = r ^ y = 

0 of particle 1. This yields a quadratic form with non¬ 
negative eigenvalues proportional to r. Thus, far out in 
the channel, the particles perform almost harmonic oscil¬ 
lations with high frequencies in the directions transverse 
to the channel while moving slowly in the direction of the 
channel. This confines the motion to a small vicinity of 
the invariant manifold. Moreover, this motion is regular. 
Close to the central region, the motion is chaotic. Hence, 
we expect our few-body system to display intermittency. 

To analyze the stability of the central region, we con¬ 
fine ourselves to computing the full phase-space mon- 
odromy matrices of certain periodic orbits within the 
manifold. All orbits up to code length three (12 in total) 
are taken into account. Together they explore a signifi¬ 
cant part of the central region. For each periodic orbit 
there are five pairs of stability exponents that correspond 
to the directions perpendicular to the invariant manifold. 
A sixth pair of stability exponents corresponds to the mo¬ 
tion within the invariant manifold. The remaining two 
pairs of stability exponents vanish because of the contin¬ 
uous symmetries of the Hamiltonian. Our computations 
show that the motion in the directions perpendicular to 
the manifold is unstable. However, the corresponding 
stability exponents are significantly smaller than the sta¬ 
bility exponents which correspond to the chaotic motion 
within the invariant manifold, see Table ^ This combina¬ 
tion of large stability exponents within the manifold and 
small stability exponents perpendicular to the manifold 
leads us to conjecture that the manifold may scar wave 
functions of the corresponding quantum system. 

Clearly this is a particular property of the invariant 
manifold we have chosen. If we consider the collinear 
configuration with mirror symmetry about the center of 
the system, we obtain a different invariant manifold with 
very large transverse Lyapunov exponents. On this in¬ 
variant manifold we do not expect scarring. 

Quantum computation. To prove our conjecture we 
consider the time evolution of a Gaussian wave packet 


^'(r, t) = cexp 



ro) A(r - ro) -k -Po (r - ro) 


( 2 ) 


where < p >= po and < r >= ro define a point on the 
invariant manifold. We have used the shorthand notation 
r = (ri,r 2 ,r 3 ,r 4 ) and p = (pi, P 2 , Pa, P 4 ) for configu¬ 
ration and momentum space vectors, respectively. The 
autocorrelation function C{t) =< '^{t = 0)|'I'(t) > is 
computed in semiclassical approximation. Within the 
manifold we used Heller’s cellular dynamics [|^ which 
takes into account the nonlinearity of the classical mo¬ 
tion. In the transverse direction the time-propagation 


was done using linearized dynamics only. This approxi¬ 
mation neglects any recurrences from the transverse di¬ 
rections and implies a permanent flux of probability out 
of the manifold and its vicinity. On the time scales con¬ 
sidered here, the linearization is justified since the clas¬ 
sical return probability to the manifold of transversely 
escaping trajectories is negligible. It is also important 
to note that the loss of probability inside the manifold 
is not severe since the transverse stability exponents are 
not too large. 

We launch wave packets along periodic or aperiodic 
orbits lying within the invariant manifold and consider 
their revival as measured by the autocorrelation function. 
To achieve shorter recurrence times the initial packet was 
symmetrized with respect to the reflection symmetry of 
the system within the invariant manifold. This effects 
a partial projection onto the symmetric and eliminates 
the antisymmetric representation of ^ 4 ; antisymmetriz¬ 
ing with respect to the reflections we could have elimi¬ 
nated the symmetric rather than the antisymmetric rep¬ 
resentation. 

We propagate such wave packets for both, the man¬ 
ifold with small and with large transversal exponents. 
For not too unstable periodic orbits we expect a fairly 
strong revival after one period, known as the linear re¬ 
vival 1^. As an example, we show in Fig. ||the real part 
of the autocorrelation function, calculated for the more 
stable invariant manifold and starting on the periodic or¬ 
bit 2. We indeed find strong linear revival. However, at 
larger times we find randomly scattered strong revivals, 
the revival corresponding to the second period not being 
dominant. This implies that a significant fraction of the 
original amplitude remains within the invariant manifold, 
and that this fact is not related to the periodic orbit we 
started on. Revivals calculated for packets started on 
aperiodic orbits show similar features except for the ob¬ 
vious absence of the linear revival. The collinear and 
more unstable of the two manifolds discussed above, on 
the contrary, shows practically no revival at all. Even the 
linear revival is negligible, see Fig. ||, despite the fact that 
the packet was launched along a short periodic orbit. 

Fig. ^ shows the Fourier transform of the autocorre¬ 
lation function whose real part was displayed in Fig. H. 
Our time sequences are not long enough to resolve the 
spectrum, and it seems difficult to obtain sequences of 
sufficient length with a reasonable numerical effort. Nev¬ 
ertheless the Fourier transform in Fig. ^ is very informa¬ 
tive. It shows a peak in energy (frequency) corresponding 
to the energy of the original wave packet. We have chosen 
the initial conditions (po, ro) such that the energy of the 
initial wavepacket is centered around E = 1. An integra¬ 
tion over classical phase space shows that this energy is 
chosen around the level. This implies that we are 

deep within the semiclassical region. The expected peak 
in energy has a superposed equidistant structure with 
spacing AE = 0.034 corresponding to the period of the 
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orbit on which we launched the wave packet. In addition 
there is a significant superposition of other frequencies 
resulting in a finer, non-equidistant structure. 

This indicates that other trajectories in the manifold 
contribute. The structure under the peak is typical for 
the more stable invariant manifold and practically absent 
for the collinear one (not shown). 

Our argument suggests that the enhanced revival seen 
on sufficiently stable invariant manifolds should not de¬ 
pend on the irreducible representation of ^4 we choose, 
unless the invariant manifold itself limits the possible rep¬ 
resentations. We could not check this point since our cal¬ 
culation did not project onto a specific irreducible rep¬ 
resentation. Nevertheless, the localization effect implicit 
in the large revival of the autocorrelation function is not 
of classical nature. Indeed, an ensemble of trajectories 
launched near the invariant manifold did not return near 
it on the time scale we consider. We note that this obser¬ 
vation also justifies the linear approximation for cellular 
dynamics in directions transversal to the invariant man¬ 
ifold. 

In summary, we have shown that wave packets may 
have unusually long life times on certain invariant mani¬ 
folds characterized by small classical transverse instabil¬ 
ities. This is a quantum and not a classical phenomenon 
and constitutes a novel and very exciting extension of 
the concept of a scar to which we attribute considerable 
significance. Indeed, invariant manifolds of the type con¬ 
sidered above occur generically in many-body systems of 
physical interest like nuclei, atoms, molecules, or small 
metallic clusters. The importance of these manifolds for 
quantum properties will hinge on their stability proper¬ 
ties: Stability will determine the degree to which scars 
actually exist in such systems. Stability is a system- 
specific property, of course, and generic statements are 
at least very difficult. We recall, however, that in Helium 
the collinear manifold (which is linearly stable in the 
perpendicular direction) supports doubly excited states. 
This shows why we expect our results to be generic for 
interacting few-body systems. 
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Code 

T 

M|| 


n 

1 

3.313 

5.74 

— 

0 

2 

2.622 

4.86 

3.29 

2 

01 

2.958 

4.46 

1.22 

2 

02 

4.036 

6.99 

4.18 

3 

12 

7.933 

11.67 

— 

0 

001 

8.015 

9.64 

— 

0 

002 

5.074 

8.07 

6.84 

4 

Oil 

4.541 

7.30 

0.35 

2 

012 

5.753 

10.00 

1.63 

1 

022 

6.623 

11.37 

7.13 

5 

112 

2.860 

4.79 

0.68 

1 

122 

13.259 

23.03 

21.41 
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TABLE I. Periods and stability exponents of periodic or¬ 
bits up to code length three lying inside the invariant man¬ 
ifold. T is the period, M|| the stability exponent inside the 
manifold, total) real positive 

stability exponents perpendicular to the manifold. The en¬ 
ergy of the orbits is E = 4. 


( 4 ) ( 1 ) 


( 3 ) ^ ( 2 ) 

FIG. 1. Collective configuration 
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FIG. 2. Autocorrelation function C{t) of a symmetrized 
wave packet launched on a periodic orbit with period T inside 
the weakly unstable manifold. In addition to the linear revival 
around t = ^, a strong nonlinear revival is seen for larger 
times. 
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FIG. 3. Autocorrelation function C{t) of a symmetrized 
wave packet launched on a short periodic orbit with period 
T inside the very unstable collinear manifold. Practically no 
revival is seen. 
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FIG. 4. Fourier transform S{E) of the autocorrelation 
function C[t) shown in Fig. ^ The equidistant structure cor¬ 
responding to half the period of the periodic orbit is accompa¬ 
nied by a fine structure resulting from the nonlinear revival. 
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